Chain End Mobilities in Polymer Melts - A Computational Study 
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The Rouse model can be regarded as the standard model to describe the dynamics of a short 
polymer chain under melt conditions. In this contribution, we explicitly check one of the fundamental 
assumptions of this model, namely that of a uniform friction coefficient for all monomers, on the 
basis of MD simulation data of a PEO melt. This question immediately arises from the fact that 
in a real polymer melt the terminal monomers have on average more intermolecular neighbors than 
the central monomers, and one would expect that exactly these details affect the precise value of the 
friction coefficient according to the Fluctuation-Dissipation Theorem. The mobilities are determined 
by our recently developed statistical method 1 , which provides complementary information about the 
local polymer dynamics and thus contrasts standard quantities such as the mean square displacement 
(MSD) or the Rouse mode analysis. Within our analysis it turned out that the Rouse assumption of 
a uniform mobility is fulfilled to a good approximation for the PEO melt, although the underlying 
microscopic dynamics is highly affected by different contributions from intra- and intermolecular 
excluded volume interactions, which, each in itself, cannot be taken into account by a modified 
friction coefficient. However, in the sum, these non-trivial terms roughly cancel each other, and 
the overall dynamics corresponds to that of a phantom chain. These effects remain elusive when 
studying the dynamics with the MSD only, since this observable characterizes a sum over dynamical 
contributions due to the bare mobility, chain connectivity and more complicated local potentials. 

PACS numbers: 36.20.Ey 61.25.H- 
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I. INTRODUCTION 

The Rouse model^ is one of the standard models to 
describe the dynamics in a polymer melt of non-entangled 
chains. Here, the polymer chain is modeled as a se- 
quence of N harmonically linked beads. All intermolec- 
ular interactions of the chain are reduced to a frictional 
and a stochastic force, both characterized by the friction 
coefficient ( via the Fluctuation-Dissipation Theorem- 
Within this model, £ takes the same value for all beads 
irrespective of the monomer position n. However, when 
going to real polymer melts, it is obvious that the indi- 
vidual polymer segments along the chain cannot have ex- 
actly the same intermolecular environment due to chain 
connectivity. A terminal segment has on average more in- 
termolecular neighbor segments in its first coordination 
sphere than a segment located in the center of the chain. 
Therefore, it is questionable if the assumption of a uni- 
form friction constant along the entire chain as assumed 
in the Rouse model still holds for a realistic chain in the 
melt. 

The dynamics of the chain ends also plays a special role 
in the limit of long chains, which is usually described by 
the reptation model^. Within this picture, the topologi- 
cal constraints imposed by the other chains are modeled 
as a tube-, and the tagged chain performs Rouse-like 
motion within this effective tube. Here, the local dy- 
namics of the chain ends play an important role in a 
twofold manner: First, since the chain can only escape 
the tube at its ends, one would expect that the dynamics 
of the chain ends significantly influences the overall re- 



laxation mechanism. For example, it has been observed 
in simulations^ — that the tube constraints are less pro- 
nounced for the terminal monomers. Improved models 
suggest that the tube relaxation is enhanced by so-called 
contour-length fluctuations^, which have also been ob- 
served experimentally^. Here, the motion of the chain 
ends leads to a loss of memory of the initial tube. Sec- 
ond, switching to the surrounding chains imposing the 
tube constraint, it was argued that the chain ends do 
not contribute to the formation of entanglements, and a 
revised tube model has been proposed*^. Experimen- 
tal data of a bimodal melt was successfully interpreted 
within this concept^. 

A different but related interplay of topological con- 
straints and the chain end dynamics becomes important 
if a polymer melt approaches the glass transition tem- 
perature T g . The dependence of T g on the molecular 
weight is usually described by the empirical Flory-Fox 
equation- 2 . Within this picture, the chain ends expe- 
rience a larger amount of free volume than the central 
monomers, leading to a lower T g for short chains, which 
has also been observed in MD simulations^. 

For all these reasons, the specific role of the chain-end 
dynamics as compared to the motion of the central parts 
of the chain is a long-standing issue in polymer science, 
which has naturally already been investigated by several 
experimental^— and numerical^— ^ studies. However, 
when discussing 'segmental mobilities', one should keep 
in mind that this expression may refer to essentially two 
different microscopic quantities: On one hand, the seg- 
mental mobility may be equal to the inverse Rouse fric- 
tion coefficient, C _1 j which is determined by the inter- 
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molecular collisions on short time scales according to the where the R n correspond to the position vectors of the 

Fluctuation-Dissipation Theorem^. On the other hand, beads. Here, the value of kg/(k^T) = 2.1 has been cho- 

the same term is used to describe how fast a given seg- sen, for which we found that the characteristic ratio of the 

ment moves, and thus rather corresponds to a supcrpo- SFCM matches that of the PEO chains in the MD simu- 

sition of chain-connectivity effects (which are naturally lations. Naturally, the semiflexible chains contained the 

less pronounced for terminal monomers) and the bare same number of monomers (N = 48) as the PEO chains. 

Rouse mobility Within this context, one should also For the friction coefficient (, the temperature fceT and 

keep in mind that most experiments such as spin-labeling the mean squared bond length 6 2 , unit values were used. 

technique s 14 ^ 6 ' 17 or neutron reflectivity-i^/neutron spin In this case, a time-step of At — 0.002 turned out to be 

echo&ii measurements as well as the standard observ- sufficiently small in order to propagate the system with 

ables calculated from simulation data such as the mean a simple Euler integrator, 
square displacement (MSD) ^ quantify the net move- 
ment of a polymer segment, and thus rather characterize 

the complex interplay between chain connectivity, bare III. PQ-METHOD 

mobility and more complicated many-chain effects. 

In contrast to this, we focus on the first definition given To briefly review the pq-method 1 , we start with a 
above, namely the bare Rouse mobility ( : . In particu- stroboscope view on the local Langevin dynamics of 
lar, we check the classical Rouse assumption of a uniform monomer n of a polymer chain 
C- value for all polymer segments independent of their po- 
sition within the chain (i. e. at the end or in the center) p n (t, At) = -A n (At)q n (t) + noise , (2) 
for the case of a realistic PEO chain under melt con- 
ditions. This is a highly pertinent question, since it is where the definitions 
obvious that the intermolecular environment, which de- 
termines the precise value of C within the picture of Brow- p«(t, At) = (1/At) [R„(t + At) - R„(t)] (3) 
nian motion^, is significantly different for a terminal and 
a central PEO monomer. In order to address this issue, an d 
we apply our recently developed method to determine 

segmental mobilities (pq-method^) to atomistic MD sim- 1n(t) — 2R„(t) — Rn+An(*) _ Rn-A«(*) (4) 

ulations of a PEO melt. It has been devised to extract , , „ , ._. 



this specific informaton about the local friction. The re- 
sults from the MD simulations are interpreted within the 
simpler semiflexible chain model (SFCM 19 ), in which a 
Rouse chain is augmented by an additional angle poten- 
tial, thus incorporating chain stiffness. From a technical 
point of view, additional emphasis is put on the compar- 
ison of our method with conventional dynamical quan- 
tities such as the MSD. In the course of our analysis, 
interesting chain end effects are observed, which cannot 
be captured by a modified friction coefficient, but rather 
reflect more subtle intra- and interchain interactions. 




approximately fulfilled. The characteristic ratio in turn 
was determined from the mean squared end-to-end vector 
(R 2 ) using the identity = (B%)/[(N - where 
bg is the average squared distance between two chemical 
monomers (in total N ). For PEO (C^ = 3.2), An = 4 
was obtained. Due to the choice of kg in eq. [TJ the same 
An-value was used for the SFCM. 



Additionally, the SFCMi2 has been simulated via 
Brownian Dynamics (BD) Simulations. In this model, 
a Rouse chain is augmented by a bending potential 



For our analysis, MD simulation data of a PEO melt 
from a previous study^ has been used. The system con- 
sisted of 16 PEO chains with N = 48 monomers each. 
The simulations were performed in the NVT ensemble 
with the GROMACS simulation package^ using an effec- 
tive two-body polarizable force field 2 ^. The temperature 
had been maintained at T = 450 K by a Nose-Hoover 
thermostat. 



II. SIMULATION DETAILS 



Eq. [5] serves as a definition of a time-dependent, ef- 
fective mobility A n (At), which for a Rouse chain with 
An = 1 approaches the true ratio k/Q in the limit 
At — > 0. Here, k = ik-oT/b 2 is the entropic force con- 
stant and b 2 the average squared bond length between 
two Kuhn segments. For finite At and appropriate An 
eq. [5] maps the local dynamics of a more realistic chain 
(i. e. the SFCM or PEO) onto the Langevin dynamics of 
a Rouse bead, thus yielding an effective mobility A n (At) 
for bead n on the time scale At. From simulation data, 
A n (At) can easily be determined by a linear regression 



A, (At) 




(5) 




(1) 



where a represents the Cartesian coordinates and 
(p„ lQ (At)) = (q n , a ) = due to isotropy. 
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Figure 1: Sketch depicting various definitions used in the 
present analysis for a simplified excluded volume chain. De- 
pending on the direction in which the center of mass of the 
subchain consisting of all monomers within q n moves during 
time At, one can distinguish between head (H) and tail (T) 
monomers (see text for further explanation). 



Because of the additional backdragging forces, exerted 
by the remaining part of the polymer chain at longer 
times, A n (At) naturally decreases with increasing At. 
On the other hand, the dynamics in the femtosecond 
regime is governed by the local chemical structure and by 
ballistic effects. We found in our previous worki on the 
same simulation data that A(At), when averaged over 
all monomers except the five outermost, shows quanti- 
tative agreement for the SFCM and PEO from approx- 
imately 5 ps. Due to this agreement, we will use the 
A n (At)-cmves of the individual SFCM monomers to in- 
terpret the results for PEO. Moreover, we also observed 
in refi that A(At) for both PEO and the SFCM con- 
tains additional anharmonic contributions on time scales 
shorter than several ten to hundred picoseconds. Thus, 
for the short-time regime, the microscopic polymer struc- 
ture manifests itself in the form of additional terms of 
higher order than the linear relationship given by Eq. [3] 
In order to remove these rather delicate features, we mod- 
ified Eq. [2] in the present analysis by additionally taking 
into account a cubic term, thus allowing us to extract the 
linear, Rousean contribution separately. Of course, this 
modification does not affect our method from a concep- 
tual point of view. 



IV. STRUCTURAL PROPERTIES 

We start with the radial distribution function of the 
ether oxygens of the PEO chains. This quantity has been 
computed for both the terminal (n = 1 and n — 48) and 
the central (n — 24 and n — 25) PEO monomers (fig. 0). 
The average has been performed over all chains in the 
simulation box as well as over all initial time frames. 
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Figure 2: Radial distribution function of the terminal (n — 1 
and n = 48) and the central (n = 24 and n = 25) PEO 
monomers in the melt. The average was performed over all 
chains and all initial time frames. For clarity, the curve of the 
central monomer has been divided by two in the inset. 



Up to approximately 0.4 nm, two intramolecular peaks 
around r — 0.28 nm and r = 0.36 nm can be observed, 
corresponding to the neighbor monomer (s). The occur- 
rence of two peaks indicates the existence of two preferred 
conformations (note that the root mean squared distance 
between the oxygen atoms of two bonded monomers is 
about 0.32 nm). Naturally, these peaks are approxi- 
mately twice as high for n = 24 due to chain connec- 
tivity (see inset of fig. [2 where the curve of monomer 
24 has been divided by two). When going to distances 
of about 0.4 — 0.5 nm, one indeed observes that the ter- 
minal monomer has more intermolecular neighbors in its 
first coordination shell than the central monomer. 

Of course, these observations are not very surprising, 
since the different intermolecular structure in the vicin- 
ity of terminal and central monomers seems only logical. 
However, as pointed out above, exactly these differences 
may also alter the dynamical behavior beyond the gen- 
eral differences due to the dissimilar chain connectivity, 
which are already contained in the Rouse model. 



V. MEAN SQUARED DISPLACEMENT 

In order to get a first impression of the local polymer 
motion, we start with the MSD, which can be regarded 
as a standard tool to study the dynamics in simulations. 

For the MSD of the entire PEO chains, we natu- 
rally observe the well-known subdiffusive center-of-mass 
motion on time scales shorter than the Rouse time tr 
(not shown), which has been found in several simu- 
lation a 18 ' 23 " — , experiment o 25 ' 27 and theoretical analy- 
se o 28 ' 29 . It is generally believed that the subdiffusivity 
for realistic chains is caused by a combination of the ex- 
cluded volume interactions and chain connectivity, giv- 
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Figure 3: MSD in the center-of-mass frame of the terminal 
and the central monomers for PEO (symbols) and the SFCM 
(dashed line). The inset shows the ratio of the center-of-mass 
MSD in direction of R c and the total center-of-mass MSD for 
PEO. The dashed line indicates the ideal value for isotropic 
diffusion. 
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Figure 4: Effective mobilities Ai(At) (top) and A2i(At) (bot- 
tom) for the terminal and the central monomers of PEO and 
the SFCM. The gray dashed line indicates the plateau value 
of An(At). Here, a distinction was made in which direction 
(relative to qi or q24.) the center of mass of the subchain 
consisting of all monomers within q n moved during At. 



ing rise to correlated motion of the chains 2 ^. More re- 
cently, also viscoelastic hydrodynamic interactions have 
been discussed 2 ^ as a reason for the anomalous center-of- 
mass diffusion. Apart from this, we additionally observed 
in the present analysis that the center-of-mass diffusion 
is slightly anisotropic with respect to the orientation of 
the end-to-end vector R c in the subdiffusive short-time 
regime. This is demonstrated by the inset of fig. El which 
shows the ratio of the component parallel to R c relative 
to the total center-of-mass MSD, i.e. ([(R CBl (f + At) — 
R cm (t)) ■ R c (t)] 2 )/([R cm (t + At) - R cm (t)] 2 >. In the ab- 
sence of anisotropies one would expect an ideal ratio of 
1/3. However, the actual ratio is larger for all At shorter 
than the Rouse time (tr ps 8.6 ns for PEO), where the 
deviations increase until they reach a maximum at about 
2 ns. This demonstrates that the preferential motion is 
along the primary axis of the chain, i. e. parallel to R . In 
how far the observed anisotropy is related to the fact that 
real polymer coils are not spherical, but rather stretched 
in direction of R y 30 ' 31 might be investigated more thor- 
oughly in future work. 

In case of the MSD of the individual monomers, one 
typically focuses on either the monomer-averaged MSD 
or, when testing the validity of specific scaling laws from 
e. g. the Rouse theory or the reptation model, the MSD of 
a few central monomers only£~— for which the charac- 
teristic features are most pronounced. For the monomeric 
MSDs, one generally observes that for intermediate time 
scales (i.e. tr/N 2 < At < tr) the outer monomers 
move faster than the central segments. Naturally, this 
regime is already highly affected by the connectivity con- 
straints of the chain, which are less present for the end 
monomers. This is also confirmed in fig. [21 which shows 
the MSD of the terminal and the central PEO monomers 



in the center-of-mass frame. As stated above, the outer 
monomers are faster for all At shorter than tr, mainly 
as a result of the chain connectivity. The MSD of the re- 
spective SFCM monomers are plotted as dashed lines in 
fig. [3] For both the central and the outer monomers, one 
observes that the MSD of the SFCM is lower than the 
PEO curve. These deviations can be related to the more 
complicated local potentials for the latter system. For 
these reasons, it is difficult to judge whether the Rouse 
assumption of a uniform £ also holds for the realistic PEO 
melt on the basis of the monomeric MSDs. 



VI. APPLICATION OF THE PQ-METHOD 

Alternatively, the local dynamics of the individual 
PEO monomers can be evaluated by the pq-method, 
yielding an effective mobility free from connectivity ef- 
fects. Fig. [H displays A n (At) for the terminal (n = 1, 
top) and the central (n = 24, bottom) PEO monomers. 
The corresponding A n (Ai)-curves for the SFCM with a 
uniform £-value are also shown. In agreement with the 
Rouse model, A n (At) (i. e. the prefactor of the linear 
term in Eq. [5]) converges to a plateau value for At — > 
that is proportional to (indicated by the gray dashed 
line in fig. S]). For larger At, A(At) decreases due to the 
dynamical contribution of the more remote monomers 1 . 
For both monomer positions, one observes a good agree- 
ment of the ,4(At)-curves of PEO (circles) and the SFCM 
(lines), thus demonstrating that the £- values of the indi- 
vidual PEO segments are essentially the same. As dis- 
cussed previously^, the SFCM basically captures the lo- 
cal dynamical features of the PEO chains in the melt 
already for short time scales. 

In order to identify possible deviations from the rela- 
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Figure 5: Effective mobilities A\(At) for the terminal 
monomer of a SFCM, in which (a) the terminal and (b) the 
four outermost monomers had a different friction coefficient 
(i.e. Ct = 0.9 C, Ct = 0.98 C, Ct = 1-02 C) and Ct = 1-1 C)> as 
well as for PEO. 

tion Ci = C24, we also calculated A\(At) for a SFCM, 
in which (a) the terminal monomer and (b) the four 
outermost monomers had a different friction coefficient 
than the other monomers (i.e. Ct = 0.9 C, Ct — 0.98 C, 
Ct = 1-02 C) and Ct = 1.1 C)- The results are shown 
in fig. [5] For the extended SFCM model curves, basi- 
cally two observations can be made: First, for a given Ct, 
Ai(At) is nearly the same for the SFCM with Ct at the 
end monomer only and the SFCM in which the four out- 
ermost monomers have a different C~ va l ue - Obviously, 
the differences become larger if the chain ends are less 
mobile than the rest of the chain (i.e. for Ct = 1.1 C) • 
Second, from the plateau values of Ai(At), one observes 
that the curves for which Ct deviates by only two per- 
cent from the C-values of the other monomers are clearly 
distinguishable. Thus, our method is able to resolve dif- 
ferences in the friction coefficient that are on the order of 
a few percent. The comparison of A\ (At) of PEO with 
the respective SFCM curves indicates that the PEO mo- 
bility of the terminal segment is lower by about 1 — 2 %. 
Nevertheless, one may conclude that the Rouse assump- 
tion of a constant mobility along the entire polymer chain 
is fulfilled for the PEO melt to a good approximation. 

In principle, when only comparing the model systems, 
one could also obtain the same information from the 
short-time MSDs (not shown). However, in case of a 
realistic polymer chain such as PEO, exactly this regime 
is dominated by dynamical contributions arising from 
the complicated, microscopic potentials (fig. El , and the 
model-chain behavior emerges only on longer time scales, 
for which the MSD is mainly governed by the global 
polymer motion. In contrast to this, our analysis is 
strictly local, as A(At) is only sensitive to the relaxation 
of a given monomer within the local curvature q„. Ac- 
cordingly, A(At) remains unaltered when studying other 
chain lengths than N = 48. 

Naturally, one would expect that the underlying PEO 
dynamics is much more complicated than that of the 



rather simplistic SFCM, which can already be seen from 
the fact that the anharmonicitics are significantly larger 
for PEO (not shown). Therefore, we investigate the di- 
rectional correlations between local chain conformation 
(as expressed by q n ) and dynamics (as expressed by 
p n ) in more detail, and modified our analysis in the fol- 
lowing way: For a given displacement ARcm" (At) = 
Rcnr 5) (t + At) - R cm ~ 5) (i) of the terminal Kuhn seg- 
ment (monomers 1 — 5, center-of-mass position Rem ), 
a distinction was made whether it moved in or against 
the direction of the qi -vector. In this way, one can distin- 
guish between head monomers (AR cm " 5 - ) (At)-qi(t) > 0), 
where the terminal Kuhn segment moves in front of the 
adjacent segments, and tail monomers (AR cm ~ 5 -'(At) • 
qi(i) < 0), where this segment follows the local chain 
contour. These two contributions to Ai(At) are also 
shown in fig. [4] (upper half). The lower half of fig. [4] 
shows a similar head-tail analysis for monomer 24. Here, 
the criterion to define head and tail monomers was if the 
center of mass of the subchain defined by the monomers 
n — Art, 11+ An moved in or against the direction of the 
q24-vector (see also sketch in fig. [I}. Note that in both 
cases the contribution of the head monomers is partly 
negative, which results from their definition. Naturally, 
trivial head and tail contributions can also be observed 
for the SFCM (fig. 0|) due to the bias resulting from the 
distinction between heads and tails. Interestingly, despite 
the good agreement of the average curves, the absolute 
values of both head and tail contributions for PEO are 
larger than for the SFCM. This is a consequence of the 
additional excluded volume, since for PEO all monomers 
within one Kuhn segment move much more correlated. 
Thus, the short-time displacement of a PEO monomer 
will already be dominated by the motion of the other 
monomers in the segment, whereas the SFCM monomers 
can interpenetrate each other and thus exhibit weaker 
motional correlations. As a result, the absolute values 
of A n (At) for head and tail monomers are larger for 
PEO. Remarkably, for both the terminal and the cen- 
tral monomers, the detailed interplay of intra- and inter- 
molecular effects leads to an overall cancellation of these 
terms, and the local PEO dynamics in the melt is essen- 
tially the same as for a phantom chain. 



VII. CONCLUSION 

In this contribution, we checked the fundamental 
Rouse assumption of a uniform friction coefficient on all 
monomers for a PEO melt. The mobilities were extracted 
from MD simulations using our previously developed pq- 
method 1 , which avoids the classical mode picture and 
rather employs a Langevin-like equation to characterize 
the local polymer dynamics. In contrast to the MSD, 
this procedure leads to the cancellation of the non-trivial 
terms for PEO, which arise from the additional chemical 
potentials. In order to interpret the local PEO dynam- 
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ics, we used a seminexible phantom chain (SFCMi^) as 
a reference. 

During the course of our analysis it turned out that 
the effective mobility of both the terminal and the cen- 
tral PEO monomers is essentially the same as for the 
SFCM from about At > 10 ps on, the deviations being 
on the order of a few percent. However, this agreement 
results from the nearly quantitative cancellation of the 
more complicated interactions in the PEO melt. A more 
detailed analysis revealed that the relaxation with re- 
spect to the local chain curvature expressed by q„ can 
be decomposed in two individual contributions (i. e. head 
and tail monomers), depending on the direction of mo- 
tion of the Kuhn segment under consideration. Due to 
the correlated motion in PEO arising from the excluded 
volume interactions, both head and tail contributions are 
larger for PEO. However, these contributions approxi- 
mately cancel, and the mobility is roughly the same for 
all PEO monomers. 



As reported previously 3 ^, our findings clearly demon- 
strate that the pq-method yields complementary infor- 
mation on the dynamics of systems involving macro- 
molecules. In further work one might study in how far 
the chain end effects persist or enhance when approach- 
ing the glass transition temperature (cf. the Flory-Fox 
equation). Moreover, the pq-method is also supposed to 
yield fruitful results for confined polymer melt a 33 ' 34 or 
complex polymer architectures 3 ^. 
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